AfHne parameterization of the dark sector: costraints from WMAP5 and SDSS 
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We study a set of universe models where the dark sector is described by a perfect fluid with an 
affine equation of state P = Pq + ctp, focusing specifically on cosmological perturbations in a flat 
universe. We perform a Monte Carlo Markov Chain analysis spanning the full parameter space of 
the model using the WMAP 5 years data and the SDSS LRG4 survey. The afflne fluid can either 
play the role of a unified dark matter (UDM), accounting for both dark matter and a cosmological 
constant, or work alongside cold dark matter (CDM), as a form of dark energy. A key ingredient is 
the sound speed, that depends on the nature of the fluid and that, for any given background model, 
adds a degree of freedom to the perturbations: in the barotropic case the square of the sound speed 
is simply equal to the affine parameter a; if entropic perturbations are present the efl'ective sound 
speed has to be specified as an additional parameter. In addition to the barotropic case, we consider 
the two limiting cases of effective sound speed equal to or 1. For a = c^ = our UDM model 
is equivalent to the standard ACDM with adiabatic perturbations. Apart of a trivial subcase, all 
models considered satisfy the data constraints, with quite standard values for the usual cosmological 
parameters. In general our analysis confirms that cosmological datasets require both a coUisionless 
massive and cold component to form the potential wells that lead to structure formation, and an 
effective cosmological constant that drives the late accelerated expansion. 

PACS numbers: 98.80.-k; 98.80.Jk; 95.35.+d; 95.36.+X 



I. INTRODUCTION 

In the last few years cosmological observations have 
become increasingly accurate, allowing various models 
to be tested or even ruled out. The one that currently 
seems to satisfy most observational requirements is the 
so-called concordance or ACDM model [1, 2]: a flat uni- 
verse, with ~ 4% of baryons, ~ 20% of an unknown 
weakly-interacting heavy component (or dark matter), 
leaving the remaining ~ 76% in the form of a cosmolog- 
ical constant (or vacuum energy density) responsible for 
the late time acceleration of the universe [3, 4]. How- 
ever, despite its simplicity, this model lacks solid theo- 
retical motivations and actually seems to require ad hoc 
assumptions, both on dark matter and the cosmological 
constant. Many hypotheses have been proposed to alle- 
viate the problem of the cosmological constant and the 
related one of why the density of the two unknown com- 
ponents are of the same order of magnitude in the present 
universe (the so-called coincidence problem). For an in- 
complete list see [5-13] and the recent review on dark 
energy [14]. 

Observationally, the dark sector is degenerate: by def- 
inition, dark components can be probed only through 
gravitational effects, leaving open a wide range of pos- 
sibilities regarding their nature and possible interactions 
[15]. In this paper we investigate an effective model for 
the dark sector, based on the affine parameterisation of 
the equation of state [16]: 



where px is the pressure, px is the energy density, and 
Pq and a are constant parameters; this leads to a time 
dependent equation of state parameter 



wx = 



px 



(2) 



Px = Po + apx, 



(1) 



An interesting property of this parameterisation is that 
it results in a constant energy density term mimick- 
ing an effective cosmological constant, with ^a = 
—Po/[Pc{l + cx)], plus an evolving term that can repro- 
duce a dark matter behaviour for certain choices of the 
parameter a. This allows one to either treat the affine 
fluid as a single unified dark component, or to use it to 
model dark energy alone. 

As shown in [17], when a is negative, this description 
can be seen as the attractor solution for a quintessence 
scalar field dynamics. Alternatively, when treating per- 
turbations, a barotropic affine fluid can be interpreted 
as a k-essence scalar field (naturally describing an ef- 
fective cosmological constant plus dark matter), while a 
scalar field with sound speed c^ = 1 act as a dark en- 
ergy component. In addition, an affine fluid description 
can also be interpreted as the result of two interacting 
dark components (one of them being a cold dark mat- 
ter component), as we discussed in detail in [18]. In the 
present work we explore several different cases resulting 
from the affine fluid description. We consider two classes 
of models: one where the affine fluid describes a unified 
dark component, the other containing a cold dark matter 
component as well. For each class, we also study three 
separate subcases, identified by the value of the speed 



of sound: the barotropic case, with c^g. = a, the case 
clg = 1, and the "silent" case [19, 20] with c^g. = 0. 

To study the properties of the model, we calculate the 
evolution of scalar perturbations in the afhne fluid by 
modifying the publicly available CAMB code, and set 
constraints to the parameters of the model by performing 
a Monte Carlo Markov Chain analysis using the cosmic 
microwave background (CMB) anisotropy WMAP 5 year 
data [21] and the large-scale matter distribution derived 
from the Sloan Digital Sky Survey (SDSS) Luminous Red 
Galaxy (LRG) 4 year data [22]. 

In the next section we describe our model and the the- 
oretical framework adopted, in Sec. Ill we discuss the re- 
sults of our numerical calculations and comparison with 
observations, and in Sec. IV we draw our final conclu- 
sions. 



turbed energy-momentum tensor components are 



(7) 



where 6^^' is the density contrast for the i component, 
F'*^ is the velocity, w'*) is the equation of state parameter 
(not necessary constant) and the subscript b refers to the 
background (i.e. unperturbed) quantities. 

In the next subsections we study in detail the be- 
haviour of the afflne dark component. 



B. Background evolution 



II. AFFINE FLUID MODEL 

A. General framework 

We perform our calculations in the context of a flat, 
homogeneous and isotropic universe, whose unperturbed 
evolution is described by the Friedman equation 



H' 



8ttG 



-P 



(3) 



where p is the total energy density, sum of the densi- 
ties of all the components in the universe, each of them 
satisfying a continuity equation that, in the case of non- 
interacting components, reads 



P(€) +'iH{p(i)+P(.,)) = 0. 



(4) 



According to the specific properties of each component 
one has different scaling behaviour: for example, for pho- 
tons and baryons p^ oc a"** and pb oc o~'^, respectively. 
We will refer to the decaying in time of the energy den- 
sity as "standard" behaviour; when the energy density 
grows in time, i.e. when /9(j') -I- p(j') < (the null energy 
condition is violated) , the behaviour is called "phantom" 
[23]. 

When treating perturbations of the background line 
element, we adopt the synchronous gauge [24]. The per- 
turbed metric then reads: 

ds'^ = a{Tf{dT'^ - {Sij + hij{x, T))dx'dx^) (5) 

where r is the conformal time and \hij | ^ 1 is the metric 
perturbation. We then compute the Einstein's equations 
at first order from the metric given above and from the 
perturbed energy-momentum tensor 






(6) 



where the index i runs over the components in the uni- 
verse, photons, baryons, and dark components. The pcr- 



The basic property of the phenomenological model we 
consider is the affine form of the pressure as a function of 
the density of the dark component, Eq. (1). Even if the 
EoS parameter of the dark component is not constant, a 
simple solution for the Eq. (4) exists and it is given by 

px = PA + (pxo-pA)a"'*'+"', a^-1- (8) 

Px = pxo-3Polna, a = -1. (9) 

where pxo is the density of the dark component at the 
present time (i.e. a = 1) and pA = —Po/{l + a), with 
a and Pq free parameters of the model. This density 
evolves in time in a way that can be either standard or 
phantom, depending on the particular choice of the pa- 
rameters. A full description of the background proper- 
ties of such a dark component is given in [25] . Here we 
want to stress that, in the absence of cold dark matter, 
this component should both be able to create the grav- 
itational potential necessary to form structures at high 
redshifts, and to drive the late time acceleration of the 
universe. With respect to a fiat ACDM model, we have 
an additional degree of freedom, a, which is the square of 
the barotropic sound speed, that allows us to investigate 
the effective equation of state of the clustering part of 
the component. 

Since the perturbation equations of the dark compo- 
nent will be written in terms of its equation of state pa- 
rameter, Eq. (2), it is interesting to explicitly consider 
the time evolution of wx ■ 

We first comment on the case pxo — Pa > (Fig. 1). 
In this case, if a > —1, wx{ci) evolves from the value 
a approaching the value —1; conversely, if a < —1, it 
approaches the value a moving away from w = —1. In 
the former situation, —3(1 + a) < and the dynamical 
part of the affine component dominates at early times. 
When a < — 1, then — 3(l-|-a) > 0, so that the evolving 
dark component increases in time, i.e. it has a phantom 
behaviour, becoming dominant at late times. The slope 
of the curve obviously depends on pxoj Po and a. 

Let us now consider the case when pxo — PA < 
(Fig. 2). The behaviour is opposite to the previous case. 



with the phantom evolution appearing when a > — 1. In 
this case there is a divergence of w in the past, mak- 
ing this choice of parameters more problematic. In this 
paper we will restrict the analysis only to cases with 
Pxo - PA > 0. 
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FIG. 1: Evolution of the dark component energy density (top) 
and equation of state parameter (bottom), for two values of 
a: a — —0.01 (left) and a = —1.5 (right). In both cases, 
pxo — PA > 0: in this case, a < — 1 results in a phantom 
regime, characterized by an energy density which increases in 
time. 
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FIG. 2: Evolution of the dark component energy density (top) 
and equation of state parameter (bottom): for two values of 
a: a = —0.01 (left) and a = —1.5 (right). In both cases, 
pxo — PA < 0: in this case, a > —1 results in a phantom 
regime, characterized by an energy density which increases in 
time. 



C. Fluid perturbations 



^i^) = -(l+^i'w)(%)+ 2) + 



3H 



dp(i 



dp{i) 






1 + W(i) 



dp(i)/dp(,) „ 
— ;— K d(j), 

1 + W(,) 



(10) 



fll) 



where we defined ikVi 



(i) 



^w- 



A pure barotropic fluid with a negative EoS param- 
eter has imaginary adiabatic sound speed that causes 
a runaway growth of perturbations. Not only this has 
unpleasant consequences for structure formation, but it 
also creates an instability in the set of coupled pertur- 
bation equations (10), (11). A viable way to overcome 
this inconvenience is to allow for entropy perturbations 
in the dark component, assuming that the effective speed 
of sound, sum of the adiabatic and entropic one, is pos- 
itive or null. We follow the formalism developed in the 
context of generalised dark matter [26] , where 



2 _ Spx 
dpx 
Px 



Cx,ad 



Wx 



^X.rost 



X 



Cx.ad 



PX 



a. 



(12) 
(13) 



Here Tx is a constant parameter we will not use since 
we prefer to specify the more fundamental quantity c^^; 
Sx.rest is the density contrast in the rest frame of the 
dark component, defined as 



3X,rest 



S 



X 



,aox 
a fc2 



(14) 



The fact that, in our fluid description, the effective speed 
of sound is a free parameter not tied to the behaviour of 
equation of state parameter wx , allows us to evade the 
tight constraints to unified dark matter models pointed 
out in [27] and arising from the runaway growth of per- 
turbations. 

To perform numerical predictions for the evolution of 
perturbations, we modified the publicly available code 
CAMB^ adding a new component whose perturbations 
are described by the following equations in the syn- 
chronous gauge: 



3JC 



^x 



^{l+Wx){Ox + 77) -3-(c^(,ff-Q;)(5x,rost 

Z a 



+ 



[l + w) 



Sx, 



X 



-X.oS ,2 



J^ ' (l + w)' 



k'5 



X.rest- 



(15) 
(16) 



Einstein's equations in the synchronous gauge and in 
Fourier space give the following system of coupled equa- 



^ http://cainb.info/ 



We adopt adiabatic initial conditions for the dark 
component [28, 29]. We first investigate the constraints 
coming from the CMB anisotropy power spectrum on a 
single dark component governed by an affine equation of 
state. As we already mentioned, this can account for both 
dark matter with a non-vanishing EoS parameter and 
a cosmological constant; we label this unified model as 
aDM model. The affine component can also be employed 
as a pure dark energy component, if standard CDM is 
present. We denominate this model as aCDM model. In 
addition to comparing our CMB anisotropy predictions 
with actual data from the WMAP 5 year observations, 
we improve our results by adding the SDSS dataset in 
order to remove degeneracies among parameters. 

In the next section we discuss the results obtained for 
both aDM and aCDM models. 



III. RESULTS 

A. Methods 

We performed a full analysis of the two classes of 
models arising from an affine equation of state (i.e. the 
aDM and aCDM models) using the Monte Carlo Markov 
Chain (MCMC) approach implemented in a modified ver- 
sion of the public CosmoMC software^ [30] . We span the 
parameter space defined by the baryon density, O^/i^, the 
cold dark matter density, flch^, the current expansion 
rate of the universe, Hq, the reionization optical depth, 
T, the spectral index rig and the normalisation amplitude 
As that parametrise the primordial curvature fluctuation 
power spectrum 



P(fc)-^(fc/fco)"=. 



(17) 



This results in a galaxy power spectrum Pg{k) = 
6|9/25P(fc). The affine dark component is charac- 
terised by the two parameters Ua (defined, as usual, as 
SttGpa/SHq) and a. Its effective sound speed squared 
has been fixed to three different values, namely 0, 1 and 
a, in order to consider the three possible clustering pos- 
sibilities, namely cold dark matter-like behaviour, scalar 
field limit and barotropic fiuid. We assumed a flat uni- 
verse and set a Gaussian prior on the Hubble parame- 
ter with mean value and standard deviation consistent 
with the Hubble Space Telescope Key Project, 72 ± 8 
km/sec/Mpc [31]. 

We computed the likelihood function of the data using 
the public code provided by WMAP team'^ that includes 
both the temperature and the polarisation CMB power 
spectrum (the main effect of the latter being a tighter 
constraint on the optical depth r). 



^ http://cosmologist.info/cosmomc/ 

^ http://lainbda.gsfc.nasa.gov/product/inap/dr2/likclihood_gct.cfm 



Even if at the background level the aDM model is 
equivalent to a dark matter with non-vanishing EoS pa- 
rameter plus a cosmological constant, there are differ- 
ences at the perturbation level; moreover, the difference 
is conceptual, since the aDM model treats the dark sec- 
tor as a whole, and can even be the result of interacting 
dark components [18]. We discuss this class of models in 
Sec. HI B. The aCDM models are discussed in Sec. HI C. 
Tables I and II summarise the best fit parameter values 
for the two classes of models. 



B. aDM Models 

In this section we investigate the properties of a 
single dark component described by an affine equa- 
tion of state. The parameters of this model are 
(Clhh'^, 0, T, In 10^°As, r^A, a). We expect the model with 
sound speed c^g- = 1 to be ruled out by the current cos- 
mological datasets: a quintessence scalar field able to 
drive the late time acceleration of the universe expansion 
prevents structure formation [17]. We tested our pipeline 
in the limit of standard ACDM model, i.e. for the choice 
a — Q, obtaining results that are in excellent agreement 
with the 5 year WMAP release [32] . In the following we 
describe the results obtained for the three sub-classes of 
models we analysed. 

aDM - c^Q = a We investigated the barotropic 
model, namely the one with c^g = c^^ = a, which does 
not require any assumption concerning entropy perturba- 
tions. As we mentioned earlier, this model has an equiv- 
alent description in terms of a k-essence scalar field. Our 
findings are shown in Fig. 3. 

With this choice of the sound speed we tested the equa- 
tion of state of dark matter. Our best fit model from the 
5 year WMAP CMB data has a == (8 ± 11) x 10-"* and 
Ha = 0.76 ±0.04: we confirm that an almost pressureless 
component is the most likely one. Since we know that 
the effect of a non- vanishing sound speed is to strongly 
modify the clustering properties, we investigated the con- 
straints which the matter power spectrum data put on 
this specific model. As expected, the constraint on a 
shrinks to |a| ^ 10 ^'', in excellent agreement with what 
found in [33]. For Ha we find Qa = 0.73 ± 0.02. In 
Figs. 4 and 5 the effect of even such a tiny barotropic 
EoS parameter is shown. 

aDM - c^g — The parameter likelihoods for the 
case of c^g = are shown in Fig. 6. The main difference 
with respect to the barotropic model is a weaker con- 
straint on a, due to the presence of a vanishing effective 
sound speed that cancels the pressure term in the pertur- 
bation equations, guaranteeing the clustering properties 
of the dark component. We get a = (—1.5 ± 3) x 10"'^ 
and r^A = 0.70 ±0.09. When the matter power spectrum 
is considered, the limit on the square of the barotropic 
sound speed a shrinks to (—2 ±2) x 10~^ at la level, and 
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FIG. 3: aDM - c^g = a: Parameter likelihoods computed 
for the aDM model under the assumption of barotropic fluid, 
i.e. a fluid that fulfils the relation c^q = a. The upper panel 
is CMB alone, the lower panel is CMB combined with the 
matter power spectrum. When the matter power spectrum 
is taken into account the constraints on the equation of state 
parameter are much tighter. The other parameters are fully 
consistent with the results of 5 years WMAP release. 



f^A = 0.69 ±0.05. 



aDM - c^g — 1 For completeness, we also performed 
the analysis in the weakly clustering limit, described by 
Cgfj = 1; as expected, the model fails completely in fitting 
the observational data. A fluid with a luminal speed of 
sound prevents the clustering at scales even close to the 
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FIG. 4; Matter power spectrum dependence on a. The black 
solid line is the matter power spectrum computed for a = 0, 
i.e. for the concordance ACDM model. The dashed curve is 
for the value q = — 1 x fO~®; the dot-dashed curve is for 
Of = 1 X fO~^. The perturbation instability is clear when a 



negative EoS parameter is chosen. 
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FIG. 5: To further illustrate the point, for the barotropic 
aDM model we plot against real data the power spectra for 
values of a at 2cr from the best fit. It is clear that the data 
constrain the value of a in two ways: 1) the theoretical curve 
has to fit the overall shape of the data distribution; 2) the data 
points at smaller scales pin down the value of |a|, constraining 
it to be small enough to i) give a small enough Jeans scale \j 
for a > 0, such that enough power is produced for A > Aj, and 
ii) for a < 0, to produce an explosive growth of perturbations 
only at small enough scales, again such that above the Jeans 
length, where gravity dominates against the pressure effects, 
the spectrum is undisturbed. It is clear from the figure that 
the second effect is dominant, in that it is extremely sensitive 
to the value of a. 



horizon [17]. 




FIG. 6; aDM - cfg = 0: Parameter likelihoods for the aDM 
model with sound speed c^g = 0. The upper panel is for CMB 
alone, the lower panel is for CMB combined with the matter 
power spectrum. The barotropic sound speed squared a is 
still consistent with 0, but the constraints are weaker than 
in the case of a pure barotropic fluid. The other parameters 
do not change significantly with respect to the concordance 
model. 



C. aCDM Models 



constant from a more general dynamical field. In this 
framework the most natural value for the speed of sound 
is c^g = 1: with this choice, our fluid description repre- 
sents well the attractor dynamics of a quintessence scalar 
field, whena<0 [17]. 



CDM - c^g = 1 In Fig. 7 we show the results for 



off' 

the aCDM model with c^g = 1. Also this case, as the 
previous one, has an equivalent description in terms of a 
scalar field, but with a standard kinetic term. The main 
effect of dark energy is to modify the low multipoles re- 
gion of the CMB power spectrum, unfortunately the one 
where high cosmic variance prevents a precise determi- 
nation of the cosmological parameters. Even worse, the 
model suffers of intrinsic degeneracies. At zeroth order, 
i.e. in the background, 1) if a ~ the dynamical part of 
the affine component behaves like dark matter, while, 2) 
if a ~ —1 it can replace the cosmological constant. Since 
we fixed the speed of sound equal to 1, the first degen- 
eracy is not present because dark matter and the affine 
component are different at the perturbation level. We are 
left with the second degeneracy, that is clearly visible in 
the flat likelihood for fl^ and the broad likelihood for a. 
We get a rather loose constraint on a, i.e. a = —1.2 ±0.4, 
while r^A — 0.5 ± 0.2. When we add the matter power 
spectrum, the CIa — a degeneracy is partially removed. 
The result is a slightly tighter constraint on a, which is 
a = -1.1 ±0.2. 



aCDM - c^jj — The results are only marginally 
affected by the value of the sound speed of the dark com- 
ponent (fig. 8), since the CMB is basically insensitive to 
the sound speed of dark energy. We find a = — 1.1±0.4 
and ilA = 0.5 ± 0.2. When the matter power spectrum is 
included in the analysis these change to a = — 1.0 ± 0.3 
and Qa = 0.3 ± 0.2 

aCDM - c^ff — a When the dark component is 
forced to be barotropic the only degeneracy we are left 
with the first degeneracy mentioned above, since a ~ 0. 
The result is flA — 0.76 ± 0.03 while flc and a are 
badly constrained {a = (6 ± 9) x 10~^): a lower value 
of Qc can be balanced by the dynamical part of the 
affine component. When the matter power spectrum 
is added we obtain a slightly tighter constraint on Q\ 
(r^A = 0.74 ± 0.02), while fie is determined by the shape 
of the spectrum. This implies actually a broad likelihood 
for the parameter a, since the coefficient (pxo ~ Pa) — 0. 
We obtain a == (1.9 ± 1.4) x 10"^. The figure 9 sum- 
marises the results described above. 



In what follows we present the results we obtained for 
the aCDM model, i.e. when we consider a fiat universe 
filled with baryons, cold dark matter and a dark energy 
component described by the affine equation of state Eq. 
(1). The choice can help to distinguish the cosmological 



IV. CONCLUSIONS 

We studied the effect of an affine EoS fiuid model ap- 
plied to the dark sector, both as a unified description 
of dark matter and an effective cosmological constant. 
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FIG. 7: qCDM - c^ff = 1: Parameter likelihoods computed 
for the qCDM model when the sound speed is fixed to c^g — 
1. The upper panel is CMB alone, the lower panel is CMB 
combined with the matter power spectrum. The almost flat 
likelihood for Oa together with the broad one for a reflect the 
degeneracies of the model. Adding matter power spectrum 
data helps to break this degeneracy since it forces Qch^ to 
be of the order of 0.11 and a ~ —1. However, Oa remains 
essentially unconstrained. 



and as a pure dark energy component. Our model makes 
use of a dynamical parameterisation relating p and p, 
as opposed to the usual cynematical description of the 
EoS parameter in terms of its current value and its first 
derivative. In a previous paper [25] we carried on a com- 
parison of the background evolution of this model with 



FIG. 8: aCDM - c^g = 0: Parameter likelihoods for the 
aCDM model when the sound speed is fixed to c^g — 0. The 
upper panel is from CMB alone, the lower panel is from CMB 
combined with the matter power spectrum. The results are 
very close to those obtained in the case of sound speed equal 
to 1. 



existing cosmological observations. In the present work, 
we focused on the behaviour of cosmological perturba- 
tions, and compared the theoretical predictions with the 
CMB WMAP 5 year data, and with the SDSS large scale 
structure data. 

As a first result, we obtained much tighter constraints 
on the parameters of the model with respect to the anal- 
ysis carried on the background observables in [25], con- 
firming that perturbations should be properly included 



TABLE I: Best fit parameter values for qDM. 



models/params 


ACDM + 
SNe & BAO (WMAP5) 


aDM - bar 
CMB MPS 


aDM - 
CMB 


MPS 


Qbh' 


0.02265 ± 0.00059 


0.0223 ± 0.0007 


0.0223 ± 0.0006 


0.0224 ± 0.0006 


0.0224 ± 0.0003 


^ch' 


0.1143 ±0.0034 


- 


- 


- 


- 


Ho 


70.1 ±1.3 


75 ±5 


71 ±2 


69 ±6 


67 ±4 


T 


0.084 ±0.016 


0.090 ±0.018 


0.083 ± 0.016 


0.086 ±0.017 


0.085 ±0.017 


Us 


0.960 ±0.014 


0.99 ±0.03 


0.960 ± 0.014 


0.959 ±0.015 


0.957 ±0.014 


log(10^M,)|fc=o.002 


3.20 ±0.08 


3.10 ±0.08 


3.19 ±0.04 


3.19 ±0.05 


3.19 ±0.04 


Oa 


0.721 ±0.015 


0.76 ± 0.04 


0.73 ± 0.02 


0.70 ± 0.09 


0.69 ± 0.05 


a 


- 


(8 ±11) X 10"'' 


(0.2 ±4) X 10" ' 


(-1.5 ±3) X 10"-' 


(-2 ±2) X 10"-' 



TABLE IL Best fit parameter values for aCDM. 



models/params 


aCD 
CMB 


M - bar 

MPS 


aCDM 
CMB 


- cis = 
MPS 


aCDM 
CMB 


- cis = 1 
MPS 


Qbh' 


0.0224 ± 0.0007 


0.0220 ± 0.0007 


0.0223 ± 0.0006 


0.0224 ± 0.0006 


0.0224 ± 0.0006 


0.0224 ± 0.0006 


Qch^ 


0.07 ±0.03 


0.106 ±0.004 


0.109 ±0.006 


0.107 ±0.006 


0.109 ±0.006 


0.107 ±0.006 


Ho 


74 ±4 


72 ±2 


71 ±6 


73 ±3 


73 ±5 


74 ±3 


T 


0.087 ±0.017 


0.083 ±0.017 


0.085 ± 0.018 


0.088 ±0.017 


0.086 ±0.016 


0.085 ±0.018 


Us 


0.975 ± 0.019 


0.964 ± 0.015 


0.962 ± 0.014 


0.964 ±0.014 


0.963 ± 0.014 


0.963 ±0.015 


log(10"'^s)|fe=0.002 


3.15 ±0.06 


3.19 ±0.05 


3.18 ±0.05 


3.17 ±0.04 


3.18 ±0.05 


3.17 ±0.05 


fiA 


0.76 ±0.03 


0.74 ± 0.02 


0.5 ±0.2 


0.3 ±0.2 


0.5 ±0.2 


0.2 ±0.2 


a 


(6 ± 9) X 10"'' 


(1.9 ±1.4) X lO"'' 


-1.1 ±0.4 


-1.0 ±0.3 


-1.2 ±0.4 


-1.1 ±0.2 





0.024 0,025 




0.02 0.04 0.06 0.0 








0.02 0.04 0.06 



in the calculations when developing effective models for 
the dark sector [34] . 



In the case when the fluid is treated as a unified dark 
component, we get values of the effective cosmological 
constant Q\ ~ 0.7, essentially independent of the speed 
of sound. For the equation of state parameter a, the 
constraints vary when the fluid is treated as barotropic 
(resulting in a slightly positive a) or a vanishing speed 
of sound is assumed (resulting in a slightly negative a). 
Both cases are however compatible with a = at one 
sigma confidence level. The inclusion of the matter power 
spectrum in the analysis has generally the effect of shrink- 
ing the confidence interval on the parameters, in partic- 
ular in the barotropic case, due to the effect of a on the 
Jeans length of the perturbations [17]. 

When standard dark matter is included, the effects of 
a on the clustering process is less relevant, because the 
matter-like component of the unified fluid is forced to 
mimic the cosmological constant behaviour. This is ap- 
parent from the fact that the a best flt value moves to 
a ^ — 1 which is typical of a cosmological constant. The 
constraints in the barotropic case remain quite tight, but 
get larger when the sound speed is set to zero. We also 
considered the case with a sound of speed equal to unity, 
which describes a scalar fleld behaviour. Also in this case 
the constraints on a are rather loose. 





0.68 0.7 0.72 0.74 0.76 0.7 




FIG. 9: aCDM 



CMB alone, upper panel, and 



likelihoods for the 
. The CMB alone 



CMB combined with MPS, lower panel, 

pure barotropic aCDM model (c^g — a] 

likelihoods show the degeneracy between Qch'' and Qa, being 

a close to 1. Adding the matter power spectrum Qch^ and 

f^A are better constrained, but we lost any information on a 

(this is because we assume a flat universe). 
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